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HIGHLIGHTS 


•  We  include  double  layer  capacitance  in  a  mathematical  model  of  lithium-ion  battery. 

•  We  demonstrate  that  the  usual  model  formulation  leads  to  a  physical  error. 

•  We  develop  a  new  formulation  suitable  for  double  layer  capacitance  implementation. 

•  We  compare  the  two  model  formulations  with  double  layer  capacitance. 

•  We  determinate  conditions  for  which  significant  differences  are  obtained. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  19  July  2013 
Received  in  revised  form 
12  September  2013 
Accepted  9  November  2013 
Available  online  28  November  2013 


Keywords: 

Lithium-ion  battery 
Electrochemical  dynamic  modeling 
Mass  and  charge  transport 
Double  layer  capacitance 


A  fundamental  model  of  lithium-ion  battery  cell  based  on  Newman’s  works  is  presented  in  this  article, 
including  migration  and  diffusion  of  species  and  charges  in  electrodes  and  electrolyte.  The  particularity 
of  our  work  is  that  the  double  layer  capacitance,  which  is  commonly  neglected,  is  integrated  in  the 
model.  We  especially  emphasize  how  to  include  correctly,  in  a  physical  point  of  view,  the  double  layer 
capacitance  in  the  electrochemical  lithium-ion  cell  model.  It  is  first  demonstrated  that  the  widely  used 
mathematical  formulation  of  the  model  leads  to  a  physical  error  when  the  double  layer  capacitance  is 
included,  and  a  new  formulation  free  of  this  error  is  proposed.  Then  are  given  the  conditions  for  which 
the  physical  error  mentioned  before  may  cause  appreciable  inaccuracy  in  the  prediction.  Simulations 
have  been  carried  out  to  validate  the  present  analysis,  and  the  results  obtained  show  that  the  difference 
in  cell  voltage  can  attain  several  tens  of  millivolts. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Numerous  mathematical  models  have  formerly  been  developed 
in  order  to  describe  locally  migration  and  diffusion  phenomena 
involved  in  a  lithium-ion  cell  through  coupled  partial  differential 
equations  (PDE)  involving  space  and  time.  One  of  the  earlier  studies 
[1,2]  led  to  one-dimensional  models  widely  used  nowadays,  and 
was  recently  improved  by  some  other  studies  introducing  thermal 
considerations  [3,4].  However,  the  double  layer  phenomenon  is 
very  frequently  ignored  in  this  kind  of  model,  although  it  can  be  of 
great  importance  to  be  integrated,  especially  in  a  system  approach, 
or  for  energy  losses  calculation  when  batteries  are  coupled  with 
power  electronics  converters. 


*  Corresponding  author.  Tel.:  +33  383  59  56  54;  fax:  +33  383  59  56  53. 
E-mail  address:  Stephane.Rael@univ-lorraine.fr  (S.  Rael). 


This  article  is  focused  on  including  electric  double-layer  phe¬ 
nomenon  in  an  electrochemical  model  of  lithium-ion  battery  cell. 
This  model  is  based  on  the  1-D  mathematical  representation  of 
charge  and  mass  transport  phenomena  in  electrodes  and  electro¬ 
lyte  developed  in  Newman’s  works  1,2].  In  the  first  part,  mass  and 
charge  transport  equations,  electrochemical  kinetic  equations,  and 
boundary  conditions  of  the  lithium-ion  battery  model  are  briefly 
reminded.  Then  emphasis  is  put  on  the  definition  of  electrolyte 
potential.  In  particular,  it  is  demonstrated  that  the  usual  formula¬ 
tion  of  the  model,  which  is  of  course  fundamentally  based  on  a  local 
definition  for  this  potential,  is  however  expressed  versus  an  elec¬ 
trolyte  potential  considered  at  the  average  electrolyte  concentra¬ 
tion.  We  propose  an  alternative  and  equivalent  formulation  of  the 
model,  which  explicitly  takes  into  account  the  local  electrolyte 
concentration  to  calculate  the  electrolyte  potential.  The  second  part 
of  the  paper  explains  how  to  include  the  double  layer  capacitance 
in  mathematical  models  for  lithium-ion  cells:  we  demonstrate  here 
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that  using  the  usual  formulation  of  the  lithium-ion  cell  model,  as 
done  like  in  Ref.  [5],  leads  to  a  physical  error  in  the  calculation  of  the 
capacitive  current  contribution.  Parameter  conditions  for  which 
this  physical  error  causes  significant  inaccuracy  are  determined  by 
equation  analysis.  Comparison  between  the  usual  formulation  of 
the  lithium-ion  cell  model,  and  the  present  formulation  relying 
upon  thorough  physical  equations,  is  finally  presented. 

2.  Mathematical  ID  model  of  lithium-ion  batteries 

As  depicted  in  Fig.  1,  a  lithium-ion  cell  is  composed  of  two 
composite  electrodes,  the  electrolyte,  a  porous  separator,  and  two 
current  collectors.  The  active  material  of  the  positive  electrode  is 
usually  a  transition  metal  oxide  containing  lithium,  such  as  cobalt 
(UC0O2)  or  nickel  (LiNiC^)  oxide,  or  iron  phosphate  (LiFePCH),  or 
NMC  materials  (LiNiyMnzCoi_y_z02).  The  insertion  material  of  the 
negative  electrode  is  a  carbon  material  of  coke  type,  graphite,  or 
amorphous  carbon.  The  electrolyte,  in  which  a  lithium  salt  e.g., 
lithium  hexafluorophosphate  LiPFg  is  dissolved,  ensures  ion  trans¬ 
port  within  the  battery.  The  separator  is  an  electronic  insulator 
which  allows  to  avoid  short  circuit  between  the  electrodes.  It  must 
be  porous  to  allow  ion  transport  between  the  electrodes.  In 
lithium-ion  technology,  it  usually  contains  polyethylene  (PE)  and/ 
or  polypropylene  (PP). 

Lithium-ion  batteries  operate  according  to  the  “rocking  chair” 
principle  [6]:  during  operation,  lithium  ions  Li+  leave  the  first 
electrode  structure,  migrate  within  the  electrolyte,  and  insert  into 
the  second  electrode  structure.  The  equations  associated  with 
insertion— extraction  phenomenon  of  lithium  ions  in  electrode 
materials  are,  at  the  negative  electrode: 

discharge 

Liy+dyCe  c  LiyC6  +  dy-Lih  +dy-e  (1) 

charge 

and  at  the  positive  electrode: 

discharge 

LizM02  +  dz-Li+  +  dze-  ~  Liz+dzM02  (2) 

charge 

Coefficients  y  and  z,  with  (0  <  y  <  1)  and  z  (0  <  z  <  1),  are  the 
lithium  insertion  rates  in  the  negative  and  the  positive  electrodes, 
respectively. 

In  this  section,  it  is  first  reminded  the  mathematical  lithium-ion 
battery  model  introduced  by  Doyle  et  al.  for  metal  lithium  cells  [  1  ], 
and  later  for  dual  lithium  ion  insertion  cells  [2  .  This  model  was 
recently  improved  by  Smith  and  Wang  by  experimental  validation  in 
Ref.  [7],  and  by  thermal  modeling  in  Ref.  [3].  Then,  an  equivalent 
formulation  of  the  model  is  presented.  However  in  the  present  work, 
the  local  electrolyte  concentration  has  been  explicitly  taken  into 


account  for  calculation  of  the  electrolyte  potential,  whereas  in  the 
usual  formulation  developed  in  the  above  cited  works,  equation 
transformations  (detailed  in  Section  2.3)  leads  to  introduce  an  elec¬ 
trolyte  potential  defined  at  the  average  electrolyte  concentration. 

2  A.  Assumptions  and  unknowns 

Fig.  1  presents  the  three  domains  of  the  ID  lithium-ion  cell 
model:  the  composite  negative  electrode  containing  LiyC6  particles, 
the  separator  -  a  porous  dielectric  ensuring  electronic  insulation 
and  ion  migration  -,  and  the  composite  positive  electrode  con¬ 
taining  LizM02  particles.  The  electrode  structure  is  impregnated 
with  electrolyte.  The  model  formulation  is  based  on  the  following 
assumptions: 

1)  the  cell  temperature  is  constant  and  uniform  throughout  the 
cell, 

2)  in  the  electrolyte  phase,  mass  and  charge  transports  are 
respectively  due  to  diffusion  after  Fields  law  and  migration  after 
Ohm’s  law,  along  the  x  direction  in  Fig.  1.  Convection  is 
neglected, 

3)  in  the  solid  phase,  lithium  transport  is  caused  by  diffusion 
expressed  using  Fields  law  and  is  treated  spherical:  the  active 
material  particles  are  supposed  to  be  monodisperse  spheres. 
Particle  radius  is  considered  small  enough  compared  with  the 
thickness  of  both  electrodes,  so  that  the  lithium  concentration  at 
solid/electrolyte  interface  could  be  computed  as  a  continuous 
function  of  space  variable  x.  Electron  transport  through  migra¬ 
tion  (Ohm’s  law)  is  one-dimensional  along  the  x  direction  in 
Fig.  1. 

Four  physical  variables  have  been  accounted  for  in  the  model: 

1 )  the  ionic  concentration  ce  of  the  electrolyte, 

2)  the  lithium  concentration  cs  in  active  material  particles, 

3)  the  electrolyte  phase  potential  </>e, 

4)  the  solid  phase  potential  </>s. 

Electrolyte  phase  variables  ce  and  </>e  can  be  defined  in  the  three 
domains,  whereas  solid  phase  variables  cs  and  </>s  are  only  in  the 
electrode  domains.  As  mentioned  before,  in  the  model  formulation 
commonly  used  in  the  literature,  electrolyte  potential  </>e  is  calcu¬ 
lated  at  the  average  electrolyte  concentration,  corresponding  to  the 
OCV  measurement  conditions. 


2.2.  Usual  model  formulation 


The  usual  model  formulation  has  been  presented  in  Ref.  [8]. 
Table  1  hereafter  summarizes  transport  and  conservation  equations 

Table  1 

ID  electrochemical  model  equations. 


Unknown  Partial  differential  equations 


Boundary  conditions 


f  Mr) 

1  Mrl 
l  at 

(  Js(x) 

I  a/sM 
t  ax 

£e.a^) 

fje(x) 
I  3 m 


= 

_  1  .  a  (n  r2  .ac,(r)\ 
~  W  d?\Us  '  Tr  ) 

-  _ ff 

~  0  s,eff  ax 

=  -JLi(x) 


r  *s(o)  =  o 

\<Ps(Rs)  =0 


-  4  (De,eff  '^w"j  "I"  t1  -  LPJ 

_  ,T  a0e(x)  ^  fv\  aln(Ce(x)) 

—  °e,eff  W '  ax  °be,effW’  ax 
=  JLi(*) 


(Js(  0)  =  Ice  n/Aceii 
I  Js(<5 neg)  =  0 
I  Js(L cell  —  <5pos)  =  0 
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f  &(0)  =  0 
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Fig.  1.  1  -D  model  of  the  lithium-ion  cell. 
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of  the  1-D  electrochemical  lithium-ion  cell  model.  In  these  equa¬ 
tions,  jLl  is  the  volume  rate  of  Li+  current  generation  (in  A  m-3), 
determined  by  Butler- Volmer  kinetics  law  as  follows: 

jLi(x)  =  as  jo •  (exp -  exp(^  -  ar (3) 

with  jo  the  exchange  current  density,  as  the  electrode  active  surface 
per  electrode  volume  unit,  p  the  local  interface  overvoltage,  Uj  the 
thermal  voltage  (uj  =  RT/F ;  R  being  the  gas  constant,  T  the  tem¬ 
perature,  and  F  the  Faraday’s  constant),  a0  and  aT  the  oxidation  and 
reduction  transfer  coefficients,  respectively.  Electrode  particles 
being  supposed  spherical,  ratio  as  can  be  expressed  versus  particle 
radius,  Rs ,  and  active  material  volume  fraction,  es,  by: 


The  effective  electronic  conductivity  crs  eff  depends  on  the  elec¬ 
trode  conductivity,  gs,  and  the  active  material  volume  fraction  es: 


0's, eff  —  es'Cs  (5) 

The  effective  salt  diffusion  coefficient  De,e ff,  and  the  effective 
ionic  conductivity  <re,eff  of  the  electrolyte  are  given  by  Bruggeman’s 
relation  [9]: 

(Deeff  =  fDe„  (6) 

l  ^e,eff  =  £e'&e 

where  (3  is  the  Bruggeman’s  exponent,  De  is  the  salt  diffusion  co¬ 
efficient  in  the  electrolyte  phase,  and  oe  is  electrolyte  ionic  con¬ 
ductivity.  For  this  last  parameter,  which  depends  on  the  electrolyte 
concentration  ce,  the  empirical  law  proposed  in  Ref.  [3]  has  been 
used: 


<re(x)  =  0.00158  ce(x)exp(  -  5.363  10_5-c]  4(x)).  (7) 

The  model  formulation  commonly  used  in  the  literature  in¬ 
volves  an  electrolyte  potential  0e  defined  at  the  average  electrolyte 
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negative  electrode 


concentration.  As  a  result,  the  local  electrode  overvoltage  involved 
in  the  jLl  expression  (Eq.  (3))  is  written  as: 

m  =  (0s(x)-0e(x))-[/(x),  (8) 

where  the  interfacial  equilibrium  voltage  U  only  depends  on  the 
lithium  insertion  rate  at  the  surface  of  solid  phase  particles,  i.e.,  on 
the  ratio  cs>e(x)/csmax  between  lithium  concentration  at  r  =  Rs  and 
maximum  lithium  concentration  in  solid  phase.  The  variations  of 
potential  U  with  the  lithium  insertion  rate  plotted  in  Fig.  2,  has  been 
determined  by  means  of  empirical  formula  given  in  Ref.  [3].  As  a 
matter  of  fact,  as  demonstrated  in  the  next  section,  in  the  usual 
model  formulation  the  theoretical  effect  of  the  local  electrolyte 
concentration  on  electrochemical  kinetics  is  implicitly  included  in 
0e  definition,  and  in  the  expression  of  the  effective  diffusional  ionic 
conductivity  crDe,eff»  which  is  given  in  this  case  by: 

^De,eff  =  2 -Uj  •  (tp  —  1)  •<7e  eff  •  (9) 


2.3.  Model  formulation  with  explicit  electrolyte  concentration  effect 

In  the  following  formulation,  the  actual  electrolyte  potential  </>e, 
defined  at  the  local  electrolyte  concentration,  is  considered.  The 
total  ion  current  density  J  e  associated  with  ion  transport  in  the 
electrolyte  phase  is  the  sum  of  two  conduction  components  for 
potential  gradient-induced  migration  of  anion  and  cation  and  two 
diffusion  terms  for  anion  and  cation  transport: 

f  J  e  =  J  n  +  J  p 

{  fn  =  -<Tn  grad(<^e)  +  (Tn-UT-grad(ln(Ce)))  0°) 

l/p  =  -ffpgrad((/>e) -ffpUTgrad(ln(ce)) 

where  Xn  and  7p  are  the  anion  and  cation  current  densities,  and  on 
and  g p  are  the  anion  and  cation  conductivities  of  the  electrolyte, 
respectively.  The  overall  set  of  equations  results  in: 

Te  =  -ffe-grad (0e)  +  (ffn  -  <rp)  •uT-grad(ln(ce)),  (11) 

where  ae  is  the  electrolyte  ionic  conductivity: 

Gq  =  (7n  +  (7p.  (12) 

It  should  be  noticed  here  that  conduction  components  of  j^  are 
added,  while  diffusion  components  are  subtracted.  In  particular,  if 
anion  and  cation  diffusion  coefficients  are  equal  (i.e.,  if  Gn  =  op, 
taking  into  account  Einstein’s  relation),  the  total  diffusion  compo¬ 
nent  is  zero,  and  the  total  ion  current  density  in  the  electrolyte  is 
purely  ohmic.  Introducing  the  cation  transference  number  tp  =  gpI 
<re,  the  total  ion  current  density  in  the  electrolyte  can  be  expressed 
versus  potential  </>e  and  concentration  ce  as  in  Table  1 : 

fe  =  -<Te-grad(0e)-ffDe-grad(ln(ce)),  (13) 

but  with  a  different  expression  for  the  diffusional  ionic  conductivity 
than  in  the  case  of  the  formulation  commonly  used  in  the  literature 
(Eq.  (9)): 

<^De  =  2-Ur(Vp— (14) 

Through  this  relation,  which  can  be  found  for  instance  in  Ref. 
[10]  in  the  presentation  of  Planck-EIenderson  equation,  it  can 
easily  be  deduced  the  following  physical  evidence.  In  a  1:1  binary 
ionic  solution,  for  which  the  diffusion  coefficients  of  the  two  ions 
are  assumed  equal  -  corresponding  to  tp  =  1/2  -  local 
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electroneutrality  implies  that  a  concentration  gradient  produces 
identical  specific  fluxes,  and  therefore  opposite  diffusion  current 
densities,  resulting  in  a  global  diffusion  current  density  equal  to 
zero. 

Eq.  (9)  can  be  established  by  considering  the  Nernst  relation 
between  the  actual  electrolyte  potential  0e,  calculated  at  the  local 
electrolyte  concentration,  and  the  electrolyte  potential  0e,mean  used 
in  the  widespread  formulation,  with  reference  to  the  average 
electrolyte  concentration: 

0e  =  0e  mean  —  Uj-  In  ( -  V  (15) 

\Ce,mean/ 

As  a  matter  of  fact,  inserting  this  relation  in  Eq.  (13)  yields  the 
charge  transport  equation  that  is  usually  given  in  the  literature: 

Te  =  -ffe-grad(0emean) -2-uT-(tp  -  l)-<Te-grad(ln(ce)). 

(16) 

The  difference  in  the  model  formulation  does  not  solely  concern 
the  modified  expression  of  the  diffusional  ionic  conductivity  in  the 
expression  of  J  e  when  the  actual  electrolyte  potential  is  used,  as 
explained  below.  As  a  matter  of  fact,  the  local  electrolyte  concen¬ 
tration  has  to  be  taken  into  account  in  the  local  OCV  calculation, 
according  to  Nernst  law.  As  a  result,  the  local  electrode  overvoltage 
involved  in  the  jLl  expression  (Eq.  (3))  is  equal  to: 


2.4.  Summary 


In  this  section,  two  formulations  of  the  same  1-D  electro¬ 
chemical  model  of  lithium-ion  cell  have  been  presented,  both 
described  by  partial  differential  equations  and  boundary  conditions 
given  in  Table  1.  The  first,  introduced  by  Doyle  et  al.  [1,2,11  ,  is 
nowadays  commonly  employed,  in  particular  by  Smith  et  al. 
[3,7,12],  White  et  al.  [4,13],  Verbrugge  et  al.  14],  Bernardi  and  Gob 
[15].  This  formulation  is  expressed  versus  an  electrolyte  potential 
defined  at  the  average  electrolyte  concentration,  corresponding  to 
the  OCV  measurement  conditions.  As  a  result,  the  OCV  considered 
for  calculation  of  the  local  electrode  overvoltage  is  the  interfacial 
equilibrium  voltage  U,  the  effect  of  the  local  electrolyte  concen¬ 
tration  on  electrochemical  kinetics  being  included  in  the  definition 
of  the  effective  diffusional  ionic  conductivity.  These  two  variables 
are  thus  given  by: 


r  ocv(x)  =  u(x) 

{  °De,eff (*)  =  2-Uj-  (tp  -  1)  -ffe.effM  ' 


(19) 


The  second  formulation,  introduced  in  this  work,  is  expressed 
versus  the  actual  electrolyte  potential,  with  reference  to  local 
electrolyte  concentration.  Then,  the  OCV  considered  for  calculation 
of  the  local  electrode  overvoltage  takes  into  account  the  effect  of 
the  local  electrolyte  concentration  through  Nernst  law,  and  the 
effective  diffusional  ionic  conductivity  does  not  have  to  be 
corrected: 


V  =  Ws-te)- (u  +  Uj\n(-S^-)),  (17) 

\  \ce,mean/  J 

which  is  equivalent  to  the  usual  formulation  expression  (Eq.  (8)), 
i.e.,  with  notations  used  in  this  section: 

V  —  (0s  —  0e,mean)  ■  (18) 


Fig.  3.  Typical  current  waveform  generated  by  power  electronic  converters,  and 
related  cell  voltage  responses.  Simulation  carried  out  with  Smith  and  Wang’s  model 
parameters  [3],  except  for  the  exchange  current  density:  Jo.neg  =  3.6  A  m2  and 
Jo.pos  =  2.6  A  m2  for  “jo/10”  curve,  jo.neg  =  0.36  A  m2  and  j0iPos  =  0.26  A  m2  for  “j0/100” 
curve.  Initial  SOC:  100%,  temperature:  300  K,  double  layer  capacitance:  20  pF  cm-2, 
switching  frequency:  10  kHz. 


r0CV(x)  =  U(x)  +  uT.ln(^) 
l  ^De,effW  =  2  -  UT  •  (tp  -2),<7e,effW 


(20) 


Important:  as  demonstrated  before,  these  two  formulations  are 
mathematically  equivalent,  and  their  computation  with  Comsol 
Multiphysics®  actually  shows  no  numerical  difference. 

When  computing  the  model,  the  usual  formulation  can  be 
preferable,  as  far  as  it  avoids  the  logarithm  function  in  the  OCV 
definition,  which  is  numerically  an  advantage.  However,  a  physical 
problem  appears  when  the  double  layer  capacitance  has  to  be 
included  in  this  lithium-ion  cell  model.  The  related  capacitive 
current  density  depends  indeed  on  the  derivative  of  the  interfacial 
voltage  0s  -  0e,  which  requires  the  actual  electrolyte  potential. 


3.  Modeling  double  layer  capacitance 

Abundant  works  have  been  devoted  to  fundamentally  based 
models  but  only  a  few  also  integrate  the  double  layer  capacitance, 
which  is  neglected  in  most  cases.  However,  describing  electric 
double-layer  phenomenon  is  necessary  for  system  applications, 
for  which  short-time  effects  are  essential.  As  a  matter  of  fact,  in 
many  applications  e.g.,  uninterruptible  power  supplies,  or  board 
networks  in  transportation,  batteries  are  often  connected  to  power 
electronic  converters  as  a  link  -  or  energy  buffer  -  to  other 
sources  e.g.,  industrial  grid,  fuel  cells,  photovoltaic  generators, 
wind  turbines,  or  to  variable  loads  e.g.,  electrical  motors.  For  safety 
reasons,  these  converters  are  usually  current  controlled  by  high 
dynamic  regulators,  with  a  response  time  to  a  step  lower  than 
10  ms.  Therefore,  accurate  knowledge  of  the  battery  resistive  and 
capacitive  behaviors  at  high  frequency  i.e.,  above  a  few  hundred 
Hertz  is  required  to  calculate  with  accuracy  regulator  parameters 
and  stability  margins  of  the  system.  Moreover,  power  electronic 
converters  generate  battery  current  ripples  at  switching  frequency, 
so  that  RMS  and  mean  currents  may  significantly  differ  from  each 
other,  especially  at  low  mean  current  level.  Therefore  an  appro¬ 
priate  model  including  double  layer  capacitance  is  required  to 
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Fig.  4.  Losses  in  floating  conditions.  Simulation  carried  out  with  Smith  and  Wang’s 
model  parameters  [3]  for  a  6  Ah  battery,  except  for  exchange  current:  jo.neg  =  3.6  A  m2 
and  jo.pos  =  2.6  A  m2  for  “j0/10”  curve,  j0,neg  =  0.36  A  m2  and  j0,pOs  =  0.26  A  m2  for  “j0/ 
100”  curve.  Initial  SOC:  100%,  temperature:  300  K,  double  layer  capacitance: 
20  pF  cm-2,  switching  frequency:  10  kHz. 


determine  accurately  battery  energy  losses  in  power  electronic 
applications. 

Figs.  3  and  4  present  simulation  results  obtained  with  the 
mathematical  lithium-ion  cell  model  detailed  in  the  previous  sec¬ 
tion.  To  introduce  typical  current  waveforms  generated  by  power 
electronic  converters,  this  model  has  been  implemented  in  Saber®, 
a  simulation  software  widely  used  in  electrical  engineering,  by 
means  of  an  electrical  analogy  to  describe  migration  and  diffusion 
of  species  and  charges  within  the  cell.  The  implementation  method 


Table  2 

Smith  and  Wang’s  model  parameters  [3]. 


Parameter 

Negative 

electrode 

Separator 

Positive 

electrode 

Thickness  5  [pm] 

50 

25.4 

36.4 

Particle  radius  Rs  [pm] 

1 

- 

1 

Solide  phase  volume 
fraction  solide  £s  [-] 

0.580 

0.500 

Electrolyte  phase  volume 
fraction  ce  [-] 

0.332 

0.500 

0.330 

Cell  area  Aceii  [m2] 

1.0452 

1.0452 

1.0452 

Max  solid  phase  cone. 

Csmax  [mol  dm  3] 

16.1 

— 

23.9 

Insertion  rate  cs>e/csmax 
at  SOC  =  0%  [-] 

0.216 

— 

0.936 

Insertion  rate  cs>e/csmax 
at  SOC=  100%  [-] 

0.676 

— 

0.442 

Average  electrolyte 

cone.  Ce.mean  [mol  dm-3] 

1.2 

1.2 

1.2 

Exchange  current 
density  jo  [A  m2] 

36 

— 

26 

Oxydation  transfer 
coefficient  aa  [-] 

0.5 

— 

0.5 

Reduction  transfer 
coefficient  av  [-] 

0.5 

— 

0.5 

Solid  phase  Li  diffusion 
coeff.  Ds  [m2  s-1] 

2  x  lO"16 

— 

3.7  x  lO"16 

Electronic  conductivity 
cr s  [S  m_1] 

100 

— 

10 

Electrolyte  phase  diffusion 
coeff.  De  [m2  s”1] 

2.6  x  10“10 

2.6  x  10"10 

2.6  x  lO"10 

Electrolyte  ionic  conductivity 
u e  [S  m-1] 

Eq.  (7) 

Eq.  (7) 

Eq.  (7) 

Bruggeman’s  exponent  [-] 

1.5 

1.5 

1.5 

Li+  transference  number  tp  [-] 

0.363 

0.363 

0.363 

Double  layer  capacitance 

Cdl  [F  m-2] 

0.2 

— 

0.2 

is  described  in  Ref.  [8].  Table  2  provides  the  geometry,  concentra¬ 
tion  and  physical  parameters  required  for  the  model  simulation: 
these  values  coming  from  Ref.  [3]  are  related  to  a  6  Ah  lithium-ion 
cell.  Fig.  3  shows  the  typical  triangular  current  waveform  induced 
by  power  electronic  conversion  system  with  a  switching  frequency 
of  10  kFIz,  and  the  cell  voltage  response  depending  on  whether  the 
double  layer  capacitance  is  taken  into  account  or  not.  A  parametric 
study  has  been  carried  out  on  the  exchange  current  density,  in 
order  to  emphasize  the  effect  of  double  layer  capacitance  on  cell 
voltage  response.  It  can  be  noticed  here  that  the  values  used  for  the 
exchange  current  density  in  Smith  and  Wang’s  work  [3]  are  large, 
around  one  order  of  magnitude  [16],  or  even  more  [17],  higher  than 
that  used  in  some  works.  As  expected,  it  can  be  observed  that 
including  double  layer  capacitance  decreases  the  “high  frequency” 
resistance.  Indeed,  this  resistance  is  only  due  to  charge  transport 
losses  in  electrodes  and  electrolyte,  whereas  electrochemical  losses 
also  contribute  to  define  the  electrical  behavior  at  high  frequencies, 
when  the  double  layer  capacitance  is  ignored. 

The  RMS  value  of  a  triangular  current  waveform  is  given  versus 
its  mean  value  /mean  and  its  current  ripple  amplitude  At  by: 


^RMS,tri  =  JmEANM/I  +  1  7  j2 - •  (21) 

V  AZ,iMEAN 

Therefore,  the  contribution  of  current  ripples  in  ohmic  losses 
increases  when  the  mean  current  decreases.  For  this  reason,  Fig.  4 
presents  simulation  results  related  to  losses  in  floating  voltage 
conditions,  i.e.,  with  zero  mean  current.  In  the  example  treated  here 
and  with  Smith  and  Wang’s  model  parameters  [3],  losses  are  about 
1.5  times  higher  if  the  double  layer  capacitance  is  ignored.  This 
figure  depends  of  course  on  activation  losses,  thus  on  exchange 
current  densities.  With  exchange  current  densities  lower  of  one 
order  of  magnitude,  the  ratio  becomes  about  3,  and  nearly  15  with 
exchange  current  densities  lower  of  two  orders  of  magnitude. 

This  section  is  aimed  at  explaining  how  to  include  the  double 
layer  capacitance  in  the  mathematical  lithium-ion  cell  model 
detailed  in  Section  2,  as  done  in  Ref.  [5].  Then,  we  emphasize  the 
theoretical  physical  importance  of  the  electrolyte  potential 
definition. 


3.1.  Including  double  layer  capacitance  in  lithium-ion  mathematical 
model  [5] 


In  a  rigorous  manner,  the  current  density  at  the  electrolyte/solid 
active  material  interface  has  faradaic  and  non-faradaic  current 
contributions.  The  second  contribution  is  related  to  the  transient 
change  in  charge  state  of  the  double  layer  capacitance,  and  is 
expressed  in  A  m-3  as: 


jDL(x) 


Cls  ■  Cdl  • 


d(0s(x)  -0eM) 

at 


(22) 


where  Cdl  is  the  surface  double  layer  capacitance  (F  itT2)  exhibited 
at  the  electrolyte/solid  active  material  interface.  In  positive  and 
negative  electrodes,  charge  conservation  equations  including  the 
double  layer  contribution  can  be  established,  as  shown  schemati¬ 
cally  in  Fig.  5.  This  leads  to: 


f  Js(x)  =  Js(x  +  dx)  +  ju(x)-dx  + jDL(x)  -dx 
1  )e  (X  +  dx)  =  Je  (x)  +  jLl  (x)  •  dx  +  jDL  (x)  •  dx  ’ 

what  can  be  locally  expressed  in  the  form: 

I W  -  ~(ju(x)  +fL(x)) 
l  TT  =  +(/U(x)  +JDL(x)) 
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x  x+dx  x 

J - 1 - 


z  «■»  - 

<t>s(x) 

— >  Js(x+dx) 

jDL(x)dx 

asC 

'  r  L 

dl^x  kC 

^)|jLl(x)-dx 

electrolyte  ,  ,  . 
phase 

<MX) 

— ►  Je(x+dx) 

Thus,  the  discrepancy  between  the  two  formulations  relies  on 
the  electrolyte  local  concentration.  More  precisely,  according  to  Eq. 
(27),  a  significant  error  can  be  observed  in  case  that  the  electrolyte 
concentration  attains  locally  very  low  levels. 

We  study  hereafter  the  conditions  to  obtain  appreciable  differ¬ 
ences  between  the  two  models  for  the  electric  double  layer  in 
lithium-ion  batteries.  For  this  purpose,  the  ID  mathematical  model 
of  lithium-ion  battery  described  in  Section  2  has  been  implemented 
in  a  numerical  code  for  multiphysics  modeling  (Comsol  Multi¬ 
physics®).  Finite  element  method  was  used,  each  domain  was  dis¬ 
cretized  in  40  elements.  Simulation  results  presented  in  this  section 
have  been  obtained  by  using  geometry,  concentration  and  physical 
parameters  extracted  from  Smith  and  Wang’s  work  [3]  and  detailed 
in  Table  2.  These  parameters  are  associated  with  a  lithium-ion  cell, 
the  capacity  of  which  is  6  Ah. 

4.1.  Differences  between  the  two  model  formulations  for  short  times 


Fig.  5.  Schematic  representation  of  local  charge  conservation  in  both  electrodes, 
including  the  double  layer  current  contribution. 


3.2.  About  lithium-ion  cell  model  formulations 


In  Eq.  (22),  the  electrolyte  potential  is  necessarily  defined  with 
reference  to  local  electrolyte  concentration,  and  not  to  mean  con¬ 
centration.  As  a  result,  using  the  widespread  formulation  of 
lithium-ion  cell  model  when  implementing  the  double  layer 
capacitance,  as  in  Ref.  [5],  is  not  physically  totally  consistent,  and 
can  theoretically  lead  to  small  inaccuracy.  In  effect,  when  relation 
(15)  between  </>e  and  </>e,mean  (cf.  Section  2.3)  is  introduced  in  Eq. 
(22),  the  double  layer  current  contribution  can  be  written  as: 


jDL{X) 


r  9(<^S(X)  -  0e  ,mean  (X))  Uj  0Ce(x) 

Qs'Cdl  at - +  fls’CDLWat_‘ 

(25) 


Only  the  first  term  is  taken  into  account  in  the  usual  formula¬ 
tion.  It  will  be  shown  in  the  next  section  that  the  second  term  is 
often  negligible,  so  that  the  two  formulations  give  very  similar 
results.  However,  for  specific  conditions  that  lead  to  a  drastic 
decrease  in  the  electrolyte  concentration  near  the  cell  current 
collectors  (x  =  0  or  x  =  Lceii,  according  to  Fig.  1 ),  a  significant  dif¬ 
ference  between  the  two  formulations  may  appear,  resulting  in  cell 
voltage  differences  of  up  to  50-100  mV. 


4.  Focusing  on  the  differences  between  the  two  formulations 

In  what  follows,  the  subscript  “mean”  is  employed  to  identify 
variables  specifically  associated  with  the  commonly  used  model 
formulation,  expressed  versus  an  electrolyte  potential  defined  at 
the  average  electrolyte  concentration.  No  particular  subscript  is 
employed  for  the  formulation  introduced  in  this  work,  expressed 
versus  the  actual  electrolyte  potential,  calculated  at  the  local  elec¬ 
trolyte  concentration.  As  mentioned  before,  capacitive  current 
contributions  can  be  expressed  in  A  m-3  as: 

j  jDL(x)  =  as-CDL-«(MizM0) 

|  i’DL  _  n  C  ^  (</*s  (x)  — 0e,mean  (*) ) 

UmeanW  -  as*LDL  at 

the  difference  between  these  two  contributions  being: 

JmeanW  -j°L(X)  =  (27) 


The  electrolyte  concentration  is  a  state  variable  that  cannot  vary 
quickly  -  in  comparison  with  electrical  variables  -  therefore  it 
cannot  be  considered  to  decrease  to  zero  within  a  few  milliseconds. 
Thus,  one  can  predict  that  there  is  nearly  no  difference  in  short 
times  between  the  two  capacitive  currents  of  Eq.  (26):  this  means 
that  the  two  lithium-ion  battery  model  formulations  are  almost 
equivalent  in  short  times,  whatever  the  conditions  of  use. 

For  illustration,  Fig.  6  presents  the  relative  difference  obtained  on 
capacitive  currents  at  the  positive  electrode  terminal 
(i.e.,  (/DL  -  Jmean)/jDLat  x  =  Lcell)  during  the  first  10  ms  of  a  high  rate 
discharge,  here  20  C.  This  regime  has  been  simulated  with  different 
values  for  the  double  layer  capacitance  CDl,  and  for  the  electronic 
conductivity  as,p0s  of  the  positive  electrode.  It  has  been  checked  that 
the  electrolyte  phase  diffusion  coefficient  De  had  no  influence  in  short 
times.  The  capacitive  current  gap  does  not  exceed  1%  at  t  =  10  ms, 
which  is  not  significant.  Consequently,  the  resulting  difference  in  cell 
voltage  at  t  =  10  ms  is  calculated  to  be  less  than  10  pV. 

Therefore,  we  can  conclude  that  for  the  description  of  “high 
frequency”  effects,  the  two  formulations  are  almost  identical.  This 
is  not  always  the  case  in  long  times,  over  a  couple  of  seconds. 

4.2.  Differences  between  the  two  model  formulations  for  long  times 

For  long  times  and  particular  conditions,  electrolyte  concen¬ 
tration  may  decrease  to  very  low  values.  Such  situation  can  occur 


Fig.  6.  Simulation  of  a  20  C  discharge:  capacitive  current  gap  in  short  times  at  positive 
electrode  terminal  (initial  SOC:  100%,  temperature:  300  K,  curve  in  black:  Smith  and 
Wang’s  parameters  [3],  others  curves:  Smith  and  Wang’s  parameters  except  for  given 
parameter  values). 


376 


N.  Legrand  et  al.  /  Journal  of  Power  Sources  251  (2014)  370-378 


when  two  conditions  are  met:  (i)  high  rate  regime,  and  (ii)  low 
electrolyte  phase  diffusion  coefficient. 

Profiles  of  electrolyte  concentration  first  depend  on  the  current 
level  and  the  working  regime  i.e.,  charge  or  discharge.  For  a  con¬ 
stant  current  discharge,  the  electrolyte  concentration  will  be 
greater  than  the  mean  concentration  ce,mean  in  the  negative  elec¬ 
trode  where  lithium  ions  are  extracted,  and  lower  than  ce,mean  in 
the  positive  electrode  where  lithium  ions  are  inserted.  Obviously 
the  reverse  situation  is  encountered  at  constant  current  charge.  At 
the  beginning  of  the  discharge,  if  the  electrode  effective  electronic 
conductivity  crs>eff  is  higher  than  the  electrolyte  effective  conduc¬ 
tivity  a 0ieff,  the  electrochemical  reaction  takes  place  preferentially 
near  the  separator,  at  x  =  dneg  and  x  =  Lceii  -  <5p0s,  according  to  Fig.  1. 
This  property  has  been  demonstrated  in  Refs.  [18]  and  [1  ].  Thus,  in 
the  case  treated  here,  one  can  first  observe  electrolyte  concentra¬ 
tion  optima  (maximum  for  the  negative  electrode,  minimum  for 
the  positive  electrode)  near  the  separator.  During  the  discharge,  the 
location  of  these  optima  moves  slowly  toward  the  cell  current 
collectors  (x  =  0  and  x  =  Lceii,  according  to  Fig.  1 ),  where  ion  flows 
are  zero.  Fig  7  gives  as  an  illustrating  example  the  time  variation  of 
the  electrolyte  concentration  profiles  within  the  lithium-ion  cell, 
for  a  discharge  at  120  A  (20  C  rate).  It  can  be  also  observed  that,  as 
expected,  a  high  rate  regime  does  not  suffice  to  obtain  low  values 
for  electrolyte  concentration.  As  a  result,  the  cell  voltage  gap  is 
lower  than  1  mV  throughout  the  20  C  discharge. 

A  parameter  which  greatly  influences  the  electrolyte  concen¬ 
tration  profile  is  the  electrolyte  phase  diffusion  coefficient  De.  Fig.  8 
presents  for  example  the  same  simulation  as  previously,  but  carried 
out  with  De/10  (in  comparison  with  Smith  and  Wang’s  value  [3]). 
With  such  a  value,  at  a  high  rate  discharge  regime,  electrolyte 
concentration  decreases  down  to  zero  at  the  positive  electrode 
terminal.  We  obtain  here  conditions,  i.e.,  high  current  rate  and  low 
electrolyte  phase  diffusion  coefficient,  for  which  the  two  formula¬ 
tions  of  lithium-ion  battery  model  including  double  layer  capaci¬ 
tance  can  exhibit  significant  differences. 

Indeed,  for  such  conditions  and  as  shown  in  Fig.  9,  the  additional 
term  in  capacitive  current  jDL  (cf.  Eq.  (25)  or  Eq.  (27))  can  represent 
an  important  part  of  jDL,  and  can  even  be  greater  in  absolute  value 
than  jDL.  In  the  example  shown  in  Fig.  9,  this  situation  occurs  for 
t  =  19.5  s.  This  clearly  means  that  capacitive  behaviors  associated 
with  the  two  formulations  become  opposite  at  positive  electrode 
terminal,  i.e.,  at  x  =  Lceii-  As  a  matter  of  fact  it  can  be  observed  in 
Fig.  10  that  the  capacitive  voltage  </>s(Lce n)  -  0e(4eii)  given  by  our 
model  formulation  goes  on  decreasing  during  discharge  after 
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Fig.  8.  Simulation  of  a  20  C  discharge:  electrolyte  concentration  profile  within  the  cell 
(electrolyte  phase  diffusion  coefficient:  De  =  2.6  x  10  10/10  m2  s_1,  initial  SOC:  100%, 
temperature:  300  K). 


t  =  19.5  s,  whereas  the  capacitive  voltage  </>s(Lce n)  -  0e,mean(^ceii)  of 
the  usual  model  formulation  begins  to  increase  during  discharge 
for  t  >  19.5  s,  which  has  no  physical  reality. 

This  results  in  a  cell  voltage  gap  between  the  two  formulations 
that  can  become  somewhat  significant,  at  several  tens  of  millivolts. 
Fig.  11  shows  the  evolution  of  this  cell  voltage  gap  in  the  extreme 
case  treated  here  (20  C  rate  discharge,  electrolyte  phase  diffusion 
coefficient  tenfold  lower  than  Smith  and  Wang’s  value  3]).  The  gap 
reaches  50  mV  at  the  end  of  the  simulation  occurring  at  t  =  27.1  s, 
time  for  which  the  value  zero  for  the  electrolyte  concentration 
gives  rise  to  computation  error. 

To  conclude,  it  should  be  noticed  that  the  low  values  used  for 
the  electrolyte  phase  diffusion  coefficient  to  emphasize  the  dis¬ 
crepancies  between  the  formulations  of  the  lithium-ion  cell  model 
are  not  only  theoretical,  but  also  practical.  Indeed,  such  values  can 
be  obtained  at  low  temperature,  as  demonstrated  in  an  experi¬ 
mental  study  of  Valoen  and  Reimers  [19]  on  transport  properties 
of  electrolyte  employed  in  lithium-ion  technology.  In  this  work, 
the  electrolyte  phase  diffusion  coefficient  is  divided  by  14  if 
temperature  decreases  from  333  I<  to  262  K,  and  the  value  at 
262  K,  of  about  5  x  10  11  m2  s-1,  is  in  the  order  of  magnitude  of 
the  value  used  here. 


Fig.  7.  Simulation  of  a  20  C  discharge:  electrolyte  concentration  profile  within  the  cell 
(initial  SOC:  100%,  temperature:  300  K). 


Fig.  9.  Simulation  of  a  20  C  discharge:  capacitive  current  contributions  at  positive 
electrode  terminal  (electrolyte  phase  diffusion  coefficient:  De  =  2.6  x  10_1°/10  m2  s-1, 
initial  SOC:  100%,  temperature:  300  K). 
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Fig.  10.  Simulation  of  a  20  C  discharge:  capacitive  voltage  at  positive  electrode  ter¬ 
minal  (electrolyte  phase  diffusion  coefficient:  De  =  2.6  x  10 ~lo/10  m2  s_1,  initial  SOC: 
100%,  temperature:  300  K). 


behaviors  associated  with  the  two  formulations  can  become 
locally  opposite,  with  notably  an  increasing  capacitive  voltage 
during  discharge,  when  the  usual  formulation  is  used.  The  cell 
voltage  gap  between  the  two  formulations  then  reaches  a  few  tens 
of  millivolts. 
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Nomenclature 


Latin  letters 


Fig.  11.  Simulation  of  a  20  C  discharge:  cell  voltage  gap  between  the  two  formulations 
(electrolyte  phase  diffusion  coefficient:  De  =  2.6  x  10-1°/10  m2  s”1,  initial  SOC:  100%, 
temperature:  300  I<). 


5.  Conclusion 

This  investigation  was  focused  on  implementing  the  double¬ 
layer  phenomenon  in  the  1-D  electrochemical  model  of  lithium- 
ion  battery  cell  developed  in  Newman’s  works  [1]  [2].  In  partic¬ 
ular,  we  demonstrate  that  the  widely  used  formulation  of  this 
model,  expressed  versus  an  electrolyte  potential  defined  at  the 
average  electrolyte  concentration,  leads  to  a  physical  error  when 
the  double  layer  capacitance  is  included.  We  propose  an  equivalent 
formulation  of  the  model,  expressed  versus  the  local  electrolyte 
concentration  for  electrolyte  potential  calculation,  which  enables  to 
take  into  account  the  double  layer  capacitance  without  the  physical 
error  introduced  by  the  usual  formulation. 

After  these  theoretical  developments,  we  investigate  the  con¬ 
ditions  for  which  the  physical  error  mentioned  before  may  cause 
significant  inaccuracy.  Equation  analysis  enables  to  determine  that 
these  conditions  are  met  if  the  electrolyte  concentration  decreases 
locally  down  to  very  low  values,  which  requires  the  combination 
of  a  high  rate  regime  and  a  low  electrolyte  phase  diffusion  coef¬ 
ficient.  In  such  conditions,  simulation  results  show  that  capacitive 


a:  active  surface  per  volume  unit  [m_1] 

A:  active  surface  [m2] 
c;  concentration  [mol  m-3] 

Cdl-  electric  double  layer  capacitance  [F  m-2] 

D:  diffusion  coefficient  of  species  [m2  s-1] 

F:  Faraday’s  constant,  96485  [C  mol-1] 

GLl:  volume  rate  of  Li+  generation[mol  m-3  s-1] 

jo:  exchange  current  density  [A  m-2] 

jLl:  volume  rate  of  Li+  current  generation  [A  m-3] 

jDL:  volume  rate  of  double  layer  current  generation  [A  m-3] 

J:  current  density  [A  m-2] 

/:  current  [A] 

N:  charge  or  species  molar  flow  [mol  s-1] 
r:  radial  coordinate  [m] 

R:  radius  [m] 
t:  time  [s] 

tp:  Li+  transference  number  [-] 

T:  temperature  [K] 

U:  electrode  equilibrium  potential  [V] 
uf.  thermal  voltage  [V] 

V:  voltage  [V] 

x :  cartesian  coordinate  [m] 

Greek  letters 


a:  transfer  coefficient  [-] 
Bruggeman’s  exponent  [-] 

5:  layer  thickness  [m] 
e:  volume  fraction  [  -  ] 
surface  overvoltage  [V] 
a:  electrical  conductivity  [S  m-1] 
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cp:  charge  or  species  molar  flow  density  [mol  m  2  s  1] 
0:  electrical  potential  [V] 

Superscripts  and  subscripts 


cell:  battery  cell 
DL:  double  layer 
e:  electrolyte  phase 
eff:  effective  value 


init:  initial  value 
mean:  average  value 
neg:  negative  electrode 
o :  oxydation 
pos:  positive  electrode 
r:  reduction 
s:  solid  phase 

s,e:  solid/electrolyte  interface 
sep:  separator 


